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Abstract. We discuss some aspects of the non-equilibrium relaxational dynamics which occur 
after a quench at a disordered critical point. In particular, we focus on the violation of the 
fluctuation dissipation theorem for local as well as non-local observables and on persistence 
properties. 



Although critical dynamics has been a subject of study for many years [1], it was rather 
recently recognized [2, 3] that, although simpler to study than glasses, they display interesting 
non-equilibrium features such as aging or violation of the Fluctuation Dissipation Theorem 
(FDT), commonly observed in more complex disordered or glassy phases [4]. In the same 
context, relaxational dynamics at pure critical point has been the subject of numerous analytical 
as well as numerical recent studies [5]. Interestingly, it has been proposed [3] that a non trivial 
Fluctuation Dissipation Ratio (FDR) X, originally introduced in the Mean Field approach to 
glassy systems, which generalizes the FDT to non equilibrium situations, is a new universal 
quantity associated to these critical points. As such, it has been computed using the powerful 
tools of RG, e.g. for pure 0{N) model at criticality in the vicinity of the upper critical dimension 
due = 4 and for various dynamics [5] . 

On the other hand, critical dynamics display interesting 'global persistence' properties. 
Indeed, it has been shown [6] that the probability Pc{t) that the global magnetization M has 
not changed sign in the time interval t following a quench from a random initial configuration, 
decays algebraically at large time Pc{t) ~ In this context, analytical progress is made 

possible, thanks to the property that, in the thermodynamic limit, the global order parameter 
remains Gaussian at all finite times t. Indeed, for a (i-dimensional system of linear size L, M{t) 
is the sum of L'^ random variables which are correlated only over a finite correlation length 
^{t). Thus, in the thermodynamic limit L/^{t) ^ 1, the Central Limit Theorem (CLT) asserts 
that M{t) is a Gaussian process, for which powerful tools have been developed to compute 
the persistence properties [7, 8]. Remarkably, under the additional assumption that M is a 
Markovian process, Oc can be related to the other critical exponents via the scaling relation 
9c = ^ = (Ac — d + 1 — r]/2)z~^, with z and Ac the dynamical and autocorrelation [9, 10] 
exponent respectively, and rj the (static) Fisher exponent. Nevertheless, as argued in Ref. [6], 
M is in general non Markovian and thus 6c is a new exponent associated to critical dynamics. 
For the non conserved critical dynamics of pure 0{N) model, corrections to this scaling relation 
were indeed found to occur at two-loops order [11], in rather good agreement with numerical 
simulations in dimensions d = 2, 3 [12]. 



Characterizing the effects of quenched disorder on critical dynamics is a comphcated task 
and indeed the question of how quenched randomness modifies these properties has been much 
less studied. In these notes, we address these questions on one prototype of such disordered 
ferromagnets, the randomly diluted Ising model: 



H = -J2 PiPjSiSj (1) 

(ij) 



where Sj are Ising spins on a d-dimensional hypercubic lattice and pi = I with probability p and 
with probability 1 — p. For the experimentally relevant case of dimension d = 3 [13], for which 
the specific heat exponent of the pure model is positive, the disorder is expected, according 
to Harris criterion [14], to modify the universality class of the transition. For 1 — p <^ 1, the 
large scale properties of (1) at criticality are then described by the following 0(1) model with a 
random mass term, the so-called Random Ising Model (RIM) [15]: 



^(V¥')' + ^h + ^(a:)]^' + |/ 



(2) 



where ip = (p{x) and ^(x) is a Gaussian random variable ip{x)'tp{x') = AS'^{x — x') and vq, the 
bare mass, is adjusted so that the renormalized one is zero. 



1. Non equilibrium dynamics of one time quantity : initial slip exponent. 

We first focus on one-time quantity, the global magnetization M{t) 

which already carries the signatures of a non-equilibrium situation. To observe them, we study 
numerically the time evolution of M{t) when the system is quenched from an initial configuration 
with short range correlations but a finite, however small, magnetization Mq. The initial stage 
of the dynamics is characterized by an increase of the global magnetization, described by a 
universal power law [9] 

M{t) ~ Mof (4) 

At larger times, to, critical fluctuations set in the system and cause the decrease of M{t) to 
zero as M{t) ~ t-{'^-2+»;)/(2z)_ j^-^ ^^^^ simulations the system is initially prepared in a random 
initial configuration with mean magnetization Mq = 0.01. At each time step, one site is randomly 
chosen and the move Si —Si is accepted or rejected according to Metropolis rule. One time- 
unit corresponds to such time steps. In all subsequent times we measure M{t) (3) for linear 
system sizes L = 8, 16, 32 and 64. Finally data are averaged over 8 x 10^ samples for L = 8 to 
10^ samples for L = 64. In the inset of Fig. 1, we show a plot of M(t) for p = 0.8 and different 
system sizes. One sees clearly that M{t) is increasing until a time to, which is an increasing 
function of L for the sizes considered here ^ and compatible with the scaling Iq L^, and then 
decreases to zero (although the aforementioned scaling for t ^ to is not clearly seen here). By 
computing M{t) for different values oi p = 0.499, 0.6, 0.65 and 0.8, we observe corrections to 
scaling, which are know to be strong in this model [16, 17]. Following Ref. [17, 18], we take 
them into account as M{t) = t'^' gp{t) with gp{t) = A'{p)[l + B'{p)t~\ where b = 0.23(2) has 
been determined previously [17, 18], A'(p), B'{p) being fitting parameters. 

^ In the thermodynamic hmit, one expects that it reaches the asymptotic value to ~ MJ~^^^^ +{<*-2+'7)/(2z)) jgj 
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Figure 1. Rescaled global magnetization 
M(t)/gp{t) as a function t for p = 0.49, 
0.6, 0.65 and 0.8 in the log-log scale. 
The linear system size is L = 64, the 
initial magnetization is Mq = 0.01 and the 
measured exponent 9' = 0.1. Inset: Global 
magnetization M{t) as a function of time t for 
p = 0.8 and L = 16, 32 and 64. 



As shown in Fig. 1, one obtains a reasonably good data collapse of M{t)/gp{t) vs. t for the 
different values of p. After a microscopic time scale, one observes a universal power law increase 
(4), from which we get the estimate 

0' = 0.10(2), (5) 
which is in good agreement with a previous two-loops estimate ^2-ioops ~ 0.0868 [19]. 
2. Two time quantities : Aging and Violation of FDT. 

Although non-equilibrium effects can already be observed in one time quantity such as the global 
magnetization M{t), it is also interesting to study two-times, t, functions, which we now focus 
on. 



2.1. Analytical approach in d = A — e. 

We study the relaxational dynamics of the Randomly diluted Ising Model in dimension d = 4 — e 
described by a Langevin equation: 

where {({x,t)) = and {C{x,t)({x' ,t')) = 2rjT5{x — x')6{t — t') is the thermal noise and r/ the 
friction coefficient. At initial time = 0, the system is in a random initial configuration with 
zero magnetization rriQ = distributed according to a Gaussian with short range correlations 

[ip{x, t = OMx',t = 0)], = To^d'^ix - x') (7) 

Notice that it has been shown that Tq~^ is irrelevant (in the KG sense) in the large time regime 
studied here [9]. We focus on the correlation C^^^ in Fourier space and the autocorrelation Cu^ 



Ctt^ = {'P{Q,t)ip{-q,t^)) , Ctt^ = {ip{x,t)ip{x,t^)) (8) 

and the response 'R-lt^ to a small external field f{—q,tw) as well as on the local response function 
respectively defined, for t > 



6f{-q,t^) ' """^ Sf{x,t^)' 



where — and (..) denote averages over disorder and thermal fluctuations respectively. In a 
previous work [18], using the exact renormalization group equation for the dynamical effective 
action, we have computed these functions 'R-Uw ' ^tt-w one-loop order: 

^L=5-'+^+''(^)'W(i-U) , Cl^=nq-'+^^^^'Fc{q'{t-t^),t/t^) ,(10) 

with (up to order 0{y^) : ■q = 0, z = 2 + [20] and 6 = [21], 6 being related to the 

autocorrelation exponent Ac through 6 = 1 + {d — 2 + r] — \c) / z. Notice that the scaling function 
is a function of g^(t — ty^) only (an therefore the superscript eg), which, although z ^ 2 
is in agreement at this order with local scale invariance arguments [22]. In Ref. [18], we have 
computed analytically the scaling functions F-ji^ at one-loop. And we will see later that some 
interesting information can be extracted from them (see in particular section 3). Here we only 
mention their asymptotic behaviors which may be relevant for our study. Defining v = q^{t — tyj) 
and u = t/tyj a, first interesting scaling regime corresponds to 3> 1, keeping v fixed. In that 
regime, one obtains an interesting relation 

Fc{v,u) = ^Fc,oo{v) + Oiu-"") (11) 

Fcoo{v) = AcoovF'j^{v) , ^Coo = 2 + 2^^ (12) 

The first relation (11) is expected from general RG arguments [9], and it has also been checked 
for pure 0{N) models [23]. The second one (12) is not predicted by such arguments, and it 
plays a crucial role in the computation of the FDR. We notice that it has also been obtained 
in the glassy phase of disordered elastic systems in d = 2 [24] as well as in the relaxational 
dynamics near the depinning transition [25]. Whether this relation holds for the present model 
at two-loops remains an open question. 

An other interesting asymptotic behavior, relevant when we want extract information on 
local quantity from (10), corresponds to u S> 1, keeping u fixed. In that regime, F^{v) (as well 
as Fc{v,u)) decays algebraically. In particular 

F'j^{v)^v-^ , i;»l, (13) 

which is in sharp contrast with pure critical systems, e.g. pure 0{N) models, where the decay 
is actually exponential [23]. 

Prom the response and correlation, it is interesting to compute the FDR X^^^ defined as [4]: 



such that Xl^^ = 1 at equilibrium. From the scalings obtained above (10), one has directly 
Xl^^ = Fx{q^{t — tw),t/tu)). In particular, in the large t/tw limit, keeping q^(t — tw) fixed, one 
has from (12) 

Jim(X^J- = 2 + /| + C?(e) (15) 



independently of v, i.e. of (small) wave vector q, which coincides of course with the asymptotic 
value for the q = mode obtained in Ref. [26]. 



It is also interesting to study the limiting value of the FDR for different observables [27, 28] , 

and in particular to focus on local quantities. Therefore we compute the FDR associated to the 
autocorrelation Cuw (8) and the local response Ruw (9) which can be written as [23]: 



For pure critical systems, T^lt^ decays exponentially for large q^t. And thus, in the limit t ^ 1, 
the integral over the Fourier modes q in (16) is dominated by the q = mode. One thus expects 
from this heuristic argument [23] that, 

lim lim XT=^ = lim Mm (17) 

We have, however, seen that for the dilute Ising model, the response function T^t^ decays 
algebraically for q^t 3> 1 (13) so that this heuristic argument does not hold. And therefore, this 
relation (17), if true at all, is far from trivial for the present model. To verify this, we need a 
direct computation in real space. 

Having obtained the scaling functions associated to and 7^^^ for any Fourier mode g, we 
obtain Ctt^ and Rtt-u, by Fourier transform [18]. Due to the algebraic large q behavior obtained 
previously (13), one obtains logarithmic corrections to scaling 

_ ^^4^+^4^Mt-U r -K ^^ + -4cln(t-U (t_\' 

(18) 

with T(u) = 1/(1 + u) + 0(e) and where A^^. arc non- universal amplitudes. From these 
expressions (18) one can then compute the local FDR under the form 

(x-f^r^=Tx{tit^) (19) 

where J^x{u) is a monotonic increasing function of : it interpolates between 1, in the quasi- 
equilibrium regime for u — 1, and its asymptotic value for u — ^ oo given by 

lim lim (Xfr°)"^ = lim lim (X^")"^ = 2 + J— + Oie) (21) 

which shows explicitly, at order 0{y/e) that the asymptotic FDR for both the global and the 
local magnetization are indeed in the same (17). 

2.2. Monte Carlo study : Autocorrelation functions. 

Let us next present results from our Monte Carlo simulations of the relaxational dynamics of the 
randomly diluted Ising model (1) in dimension d = 3, which, as previously (see section 1) were 
done on cubic lattices with periodic boundary conditions. The system is initially prepared 
in a random initial configuration with zero magnetization, and at each time step, the sites 
are sequentially updated according to Metropolis rule. Here we present our numerical data for 
p = 0.8 (other values of p are discussed in [18]). We compute the spin-spin auto-correlation 
function defined as 

= ;^E(^^(^raU) (22) 

i 
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Figure 2. Log-log plot of auto- 
correlation function Cu^ vs. t — for 
different waiting times. 
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Figure 3. Log-log plot of f2,^'^^'"Cu^ 
as function of t/tyj for different waiting 
times. 



Fig. 2 shows the auto-correlation function Ctt^ as a function oit — tyj for different values of the 
waiting time = 2^, 2^, 2^, 2^ and 2^ at p = 0.8. One observes a clear dependence on t^, 
which indicates a non-equilibrium dynamical regime. The scaling form obtained from the RG 
analysis (18) suggest, discarding the logarithmic correction, to plot tw'^^^^^Cu^ as a function of 
t/tyu. Taking the values rj = 0.0374 from Ref. [29] and z = 2.62 from Rcf. [17], we see in Fig. 3 
that, for p = 0.8, one obtains a good collapse of the curves for different t^,. As shown in Fig. 3, 
tw Ctt^ decays as a power law, which allows to estimate the value of the decaying exponent 
Xc/z: 

^ = 1.05 ± 0.03 (23) 

We have checked [18], taking carefully corrections to scalings, that this decaying exponent (23) 
is actually independent of the dilution factor p, which supports universality of the long-time 
non-equilibrium relaxation in this model. In addition, using z = 2.6, our numerical estimates of 
9' and Ac (23) are consistent with the scaling relation \f. = d — z6'. 

3. Persistence properties. 

Let us now focus on persistence properties of the dilute Ising model at criticality. Defining the 
global magnetization M{t) as in Eq. (3) we are interested in the disorder averaged probability 
Pc{t) that the magnetization has not changed sign in the time interval t following the quench. 

3.1. Analytical approach in d = 4 — e 

In a previous publication [30], we have shown that the Exact renormalization group equation 
allows to describe the time evolution of the magnetization by an effective Gaussian process M{t): 

dtM{t) + a{t)M{t) = - f dhl^u,M{ti) + at) (24) 

Jo 

1 /6e /■* 
=-2^53(7(^-0)' , a{t) = - dti^m (25) 

where ^(x) = (x + Aq^)~^, Aq being the UV cutoff, and C(i) is an effective disorder induced 
Gaussian noise with zero mean and correlations (C(*)C(*'))e// = 2Td{t — t') + Du' with: 

T /67 

^**' = yV53^^^*"^')"^^*+^')^ ^^^^ 



The idea is then to compute Pc{t) as the persistence probabihty of the process M{t). The 
rhs of Eq. (24) clearly indicates that this process is non-Markovian. However, as Tiw (25) as 
well as Dtfi (26) are of order 0{^/e), one can use the perturbative computation of 9c around 
a Markov process initially developed in Ref. [8] and further studied in the context of critical 
dynamics in Ref. [11]. In that purpose [7], let us introduce the normalized Gaussian process 

m{t) = M{t)/^J {M'^{t))eS- Let T = Int, then m{T) is a stationary Gaussian process and its 
persistence properties are obtained from the autocorrelation function: 

1 fEe / , X — I , — 1 
1 + - y — X log — log ■ 



(m(T)m(T^))eff = e-'^(^-^»U(e^-^») , " ^ ■ 4 y 53 I - ^ + 1 --4^2 

(27) 

with fj, = (Ac — d + 1 — r]/2)/z. Under this form (27), one can use the first order perturbation 
theory result of Ref. [11] to obtain the one-loop estimate: 

A = ec-fi= \/|^ ^ ~ ^ + C(e) = 0.06968... in d = S, (28) 

where /i is the value corresponding to a Markov process. The second term in Eq. (28) is the 
first correction due to the non-Markovian nature of the dynamics. Interestingly, this correction 
is entirely determined by the non-trivial structure of the scaling function A{x), which is directly 
obtained from Fc{v,u) (10). Notice also that, at variance with the pure 0{N) models [6, 11], 
these corrections in presence of quenched static disorder already appear at one-loop order. 

3.2. Numerical simulation in d = 3 

We now turn to the results from our Monte Carlo simulations of the relaxational dynamics 
of the randomly diluted Ising model (1) in dimension d = 3, which were performed on 
cubic lattices with periodic boundary conditions. The system is initially prepared in a random 
initial configuration with zero mean magnetization Mq = 0. Up and down spins are randomly 
distributed on the occupied sites, mimicking a high-temperature disordered configuration before 
the quench. At each time step, one site is randomly chosen and the move Sj — > —Si is accepted 
or rejected according to Metropolis rule. One time unit corresponds to such time steps. 
The exponent 9c is measured numerically for cubic lattices of linear size L = 8, 16, 32 and 64. 
After a quench to Tc from the initial random configuration each system evolves until the global 
magnetization first change sign. Pc(t) is then measured as the fraction of surviving systems at 
each time t, over a number of samples which varies from 2 X 10^ for L = 8 to 2 X 10^ for L = 64. 

In Fig. 4 we present the results of Pc{t) for p = 0.8 and for different lattice sizes. According 
to standard finite-size scaling [6], one expects the scaling form Pc{t) = t~^'' f {t / L^) , where z is 
the dynamical exponent. Keeping the rather well established value of 2; = 2.62(7) [17, 18] fixed, 
9c is varied to obtain the best data collapse. The final scaled plot is shown in the inset of Fig. 4. 
This allows for a first estimate of the exponent 9c 

9c = 0.35 ± 0.01 (29) 

We have also computed the persistence probability for systems quenched from random 

configuration with a small initial magnetization Mq = 0.001. The number of up A^^p and 
down Adown spins are thus : Aup = (H-Mo)/2Aocc and Adown = Aocc — Aup. First we randomly 
distribute the Aup up spins in the occupied sites of the lattice and then fill up the rest with down 
spins. As noticed previously [12] this protocol allows to reduce the statistical noise and thus to 
study larger system sizes (this however renders the finite size scaling analysis more subtle [12]). 




Figure 4. Persistence probability Pc{t) 
plotted in the log-log scale for p = 0.8 
and different L, with Mq = 0. Inset: 
L'^'^^P^it) vs i/L^ for different L with 
z = 2.62 and 6 = 0.35. 



Figure 5. Similar to Fig. 4, but prepared 
with a non-zero initial magnetization 
Mq = 0.001. The linear system size in 
this case is L = 100. 



In Fig. 5, we plot the persistence for system size L = 100 ( the data have been averaged over 
2 X 10^ ensembles) for p = 0.499, 0.6, 0.65 and 0.8. The study of larger system size allows to 
reduce the corrections to scaling. Indeed, although in the short time scales, the straight lines 
have slightly different slopes which depends upon p, at later times the slopes varies from 0.36(1) 
for p = 0.8 to 0.35(1) for p = 0.499 : this confirms the p-independent value of 9c obtained 
previously Eq. (29). 

In order to compare this numerical value (29) with our one-loop calculation (28) one needs an 
estimate for /x. Such an estimate is needed not only for the sake of this comparison but also to 
characterize quantitatively non Markovian effects. The argument mentioned in the introduction, 
relying on the CLT, which says that the global magnetization is, in the thermodynamic limit, a 
Gaussian variable is also valid in the presence of disorder, and we have checked it numerically. 
Therefore, a finite difference A (28) is the signature of a Non Markovian process. Because of 
the relatively big error bars on Ac and z, we propose, alternatively, to express /i in terms of 
the initial slip exponent , 9' [9], using Xc = d — z9' , as // = —9' -|- (1 — r]/2)/z. Using our 
previous estimate for the initial slip exponent 9' = 0.10(2) (5), this gives /Xnum = 0.27(3) and 
our numerical estimate 

Anum = 0.08 ± 0.04 (30) 

which is in good agreement with our previous one-loop estimate (28). We also notice that these 
deviations from a Markov process are slightly larger than for the pure case [12] . 

4. Conclusion 

In conclusion, we have studied, analytically using the Exact renormalization group equation as 
well as numerically, different aspects of non-equilibrium dynamics at a random critical point. 
Interestingly, concerning the violation of FDT, although the heuristic argument of Ref. [23] does 
not hold for the present model, we have shown explicitly that the limiting FDR for the local and 
global magnetization do coincide (21). And it would be interesting to compare this perturbative 
calculation with numerical simulations. In addition, in view of recent studies [18, 31, 32, 33, 34], 
it would be also interesting to study this FDR for the relaxational dynamics following a quench 
from a completely ordered initial condition. Finally, concerning the persistence properties, we 



have shown that the RG calculation, together with the perturbative methods of Ref. [8, 11] 
allows for a rather precise estimate of 9^. And it would be interesting to extend this kind of 
approach deep inside the glassy phase of disordered systems. 
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